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Abstract 

We provide a full one-loop calculation of the self energies and tadpoles of the neutral Higgs 
bosons of the NMSSM. In addition, we compute the two-loop 0(ata s + a>ba s ) corrections to the 
neutral Higgs boson masses in the effective potential approximation. With respect to earlier calcu- 
lations, the newly-computed corrections can account for shifts of a few GeV in the light scalar and 
pseudoscalar masses, and they can also sizeably affect the mixing between singlet and MSSM-like 
Higgs scalars. Taking these corrections into account will be crucial for a meaningful comparison 
between the MSSM and NMSSM predictions for the Higgs sector. 
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1 Introduction 



The Next-to-Minimal Supersymmetric extension of the Standard Model, or NMSSM [TJ[2j[3], provides 
an elegant solution to the fi problem of the MSSM [3], i.e. the question of how to relate the higgsino 
mass parameter [i appearing in the superpotential to the soft SUSY-breaking masses of the other 
supersymmetric particles. In the NMSSM the parameter [i arises as the vacuum expectation value 
(vev) of the scalar component of an additional chiral superfield S, singlet with respect to the SM gauge 
group and coupled to the MSSM Higgs superfields H\ and H2 through a superpotential term A SH1H2 ■ 
The scalar and pseudoscalar components of the singlet superfield mix with the MSSM Higgs fields of 
matching parity, while the fermion component (singlino) mixes with the MSSM higgsinos. The quartic 
Higgs scalar interaction controlled by the new coupling A can bring the additional benefit of increasing 
the tree-level prediction for the mass of the lightest Higgs boson, which in the MSSM is bounded from 
above by the Z-boson mass, in contradiction with the lower bound from direct searches at LEP [5]. 
This increase allows for a smaller contribution to the Higgs mass from radiative corrections involving 
the top quark and its scalar partner, the stop, thus reducing the so-called little hierarchy problem of 
the MSSM, i.e. the need for the stop masses to be substantially larger than the weak scale. Another 
interesting feature of the NMSSM is the existence of scenarios in which a Higgs scalar with SM-like 
couplings to fermions and gauge bosons decays mainly into a pair of lighter scalars or pseudoscalars, 
requiring a refinement of the strategies for Higgs searches at the Tevatron and at the LHC [3 [7] . A 
SM-like Higgs scalar decaying into a pair of pseudoscalars which are in turn too light to decay into 
b quarks might even have already been produced at the LEP, and escaped detection because of its 
non-standard decay chain [8]. 

Due to the crucial role of radiative corrections in pushing the prediction for the lightest Higgs boson 
mass above the LEP bound, an impressive theoretical effort has been devoted in the past two decades 
to the precise determination of the Higgs sector of the MSSM [9]. After the early realization |10j of 
the importance of the one-loop 0(at) correction^] involving top and stop, full one-loop computations 
of the MSSM Higgs masses have been provided I12j . leading logarithmic effects at two loops have 
been included via appropriate renormalization group equations [13] , and genuine two-loop corrections 
of 0(a t a s ) QUI51HE1QI1HS1, HI El US], 0{a h a s ) 1221121] and 0(a t a b + a 2 b ) [22J have been 

evaluated in the limit of zero external momentum. All of these corrections have been implemented 
in public computer codes |23[ [23] for the calculation of the MSSM mass spectrum. More recently, a 
nearly complete two-loop calculation of the MSSM Higgs masses, including electroweak effects and 
part of the external momentum dependence, has been performed [25] , and even the leading three-loop 
effects have been computed |26| . Finally, a vast literature [27] is available on the dominant one- and 
two-loop corrections to the MSSM Higgs masses in the presence of CP-violating phases in the soft 
SUSY-breaking parameters. 

x We define a t ,b = h\ b /(47r), where ht and hi, are the superpotential top and bottom couplings, respectively. Here and 
in the following we denote for brevity as 0(ott) the one-loop corrections to the Higgs masses that are in fact proportional 
to a t mt , i.e. Q^v\. Similar abuses of notation affect the other one- and two-loop corrections. 
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In comparison with the case of the MSSM, the computation of the radiative corrections to the 
Higgs masses in the NMSSM is not quite as advanced. The one-loop contributions from diagrams 
involving top/stop and bottom/sbottom loops have been computed [28] only in the effective potential 
approximation, i.e. neglecting the external momentum in the self energies. For what concerns the 
one-loop contributions from diagrams involving chargino, neutralino or scalar loops (the contributions 
arising from gauge-boson loops are the same as in the MSSM) only the leading logarithmic terms 
have been computed [29]. Similarly, among the two-loop contributions only the leading-logarithmic 
0(ata s ) and O(a^) terms - borrowed from the MSSM results under the simplifying assumption of 
fully degenerate SUSY masses - have been taken into account so far. All of these corrections have 
been implemented in a public computer code, NMHDECAY [30j, which computes masses, couplings and 
decay widths of the NMSSM Higgs bosons. 

It is clear that, for a proper comparison between the MSSM and NMSSM predictions and for a 
precise characterization of the scenarios of refs. [6j [8], it would be desirable to compute the masses 
and mixings in the NMSSM Higgs sector with an accuracy comparable to that of the calculations 
implemented in the public MSSM codes of ref. [23]. In this paper we take a few steps in this direction. 
First of all, we provide explicit formulae for the one-loop corrections to the mass matrices of the neutral 
CP-even and CP-odd Higgs bosons. In addition, we compute the two-loop 0(a t a s + a\ ) a s ) corrections 
in the approximation of zero external momentum, adapting to the NMSSM case the techniques (and, 
in part, the results) developed for the MSSM in refs. |18| , I19| , [22 j 131]. To fully match the accuracy of 
the MSSM codes it would also be necessary to include the two-loop O(o^ + atOf, + a^) corrections. 
These corrections, however, cannot be straightforwardly adapted from the MSSM case and require a 
dedicated calculation. We leave that, as well as a detailed re-analysis of the NMSSM parameter space 
taking into account the improvements in the Higgs mass calculation, to a future publication. 

The paper is organized as follows. After this introduction, in section[5]we describe the Higgs sector 
of the NMSSM. In section [3] we describe the diagrammatic calculation of the one-loop corrections 
to the CP-odd and CP-even Higgs mass matrices. In section [5] we describe the effective potential 
calculation of the two-loop 0{ata s + aj,a; s ) corrections. In section [5] we show numerical results for 
some representative scenarios. In section O we conclude. Finally, appendix A contains the definitions 
of the couplings that enter the calculation of the one-loop radiative corrections; appendix B provides 
the explicit expressions for the one-loop self energies and tadpoles that appear in the computation of 
the Higgs boson masses; appendix C contains the formulae for the two-loop corrections in terms of 
derivatives of the effective potential; appendix D contains the definitions of the functions entering the 
two-loop results. 
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2 The Higgs sector of the NMSSM 



In the SLHA conventions [321 [33] the NMSSM superpotential for the Higgs superfields Hi, H2, S and 
the quark and lepton superfields Q,U C , D c , L , E c reads 



W 



h e H 1 LE C + h d HiQD c + h u QH 2 U C - X SH X H 2 + ^ S 3 



(1) 



where we neglect colour and generation indices. The 5f7(2)-doublet superfields are contracted by the 
antisymmetric tensor e a fe, with €12 = 1. The corresponding terms in the soft SUSY-breaking scalar 
potential are 



Ksoft = raj^iTj-Hi + rr? H% H\H2 + m 2 s S*S + m 2 Q Q^Q + m 2 L L^L + mfj u* r ur + m z D d* R d R + m% e* R e R 



rt; 



^ 2 rr* 1 



,2 ~* ~ 



,2 J* 



,2 ~* 



/i e A e HiL e* R + /i d vlrf H X Q d* R + h u A u QH 2 u* R -XA x SHiH 2 + - A K S 3 + h.c. 



(2) 



where for the scalar components of the quark and lepton superfields we define u R = U c *, d R = D c * , 
e R = E c *, Q = {uL,di) T and L = (^l^l) T - For simplicity, we take the soft SUSY-breaking trilinear 
couplings (as well as the gaugino masses) to be all real, and we neglect inter-generational squark and 
slepton mixing. 

The neutral components of the Higgs fields can be decomposed into their vevs plus their CP-even 
and CP-odd fluctuations as 



H? = Vi + ^iSi + iPi) (i = 1,2) 



S 



1 



Vs+ V2 ^ 3 + iPz) ' 



(3) 



Using the minimization conditions of the tree-level scalar potential Vq to replace the soft SUSY- 
breaking Higgs masses vr? H , mji 2 and m| with combinations of the Higgs vevs and the trilinear 
couplings, the tree-level mass matrix for the CP-even fields, (Ai'g)^ 00 , reads 



/ 



(2A 2 - f)v lV2 -Xv s A 1 



(2A 2 - g 2 )viv 2 - Xv s Ay 2X 2 v 1 v s - Xv 2 (A s + kv s ) \ 



2A 2 v 2 v s - X vi {Ay, + kv s 



(4) 



\ 2X 2 VlV s - XV 2 (Ay + KV S ) 2X 2 V 2 V S - Xvi(Ay + KV S ) XA X m ^ + KV s (A K + AKV s ) ) 

where for brevity we define Ay = A\ + kv s and g 2 = (g 2 + g' 2 )/2, g and g' being the electroweak 
gauge couplings. The CP-even mass matrix is diagonalized by an orthogonal matrix R s , such that 



hi — Rij Sj , 



(5) 



where hi, with i = 1,2, 3, are the CP-even mass eigenstates ordered by increasing mass. The tree-level 
mass matrix for the CP-odd fields, (M. 2 pf Tee ' , reads 



Xv^Ay 



vi 

Xv s Ay 



\ 



Xv s Ay Xv2(Ay-3kv s ) 
Xv s ^Ay Xvi(Ay-3kv s ) 
\ Xv 2 (Ay-3kv s ) X Vi (Ay-3kv s ) 4Xkviv 2 + XA x ^-3kA k v s J 



(6) 



3 



The CP-odd mass matrix is in turn diagonalized by an orthogonal matrix R p , such that 



Rij Pj , 



(7) 



where a% stands for (G° , A\, A2). Here, G° is the neutral pseudo-Goldstone boson, while A\ and A2 
are the other CP-odd mass eigenstates, ordered by increasing mass. A commonly adopted procedure 
(see, e.g., the first paper in ref. [2j) consists in expressing Pi and P2 in terms of G° and the orthogonal 
combination A , then defining a new 2x2 orthogonal matrix, R pl ', which rotates (A , P3) into (Ai, A 2 ). 
However, we will follow the SLHA prescription and refrain from this simplification. 

In the neutralino sector, the singlino s mixes with the neutral components of the MSSM higgsinos 
hi and h®, which in turn mix with the neutral gauginos b and w . In the formalism of two-component 
spinors, the Lagrangian contains the mass terms 



■ib 



iw 



r,0 



hi hP 2 



1 


Mi 





-g'vi/V2 


g'v 2 /V2 







( 


-ib 









M 2 


g Vl /V2 


-gv 2 /V2 









—iw° 








g'v 2 /V2 





-Xv s 


— Xv 2 






hi 






gvt/y/2 


-gv 2 /V2 


-Xv s 





—X vi 






hi 












-Xv 2 


—A vi 


2kv s 


) 


\ 


s 


) 



(8) 

where Mi and M 2 are the soft SUSY-breaking gaugino masses. The neutralino mass matrix is diago- 
nalized by a unitary matrix N, such that 



X 



0) 



where tp° stands for (—ib, —iw°, hi, hi, s) and Xi are the neutralino mass eigenstates ordered by 
increasing mass. Under our simplifying assumptions all the entries of the neutralino mass matrix in 
eq. © can be taken as real. In this case N becomes a real matrix, as long as the neutralino mass 
terms are allowed to take on negative signs. 

Finally, the charged-Higgs and chargino sectors of the NMSSM are not directly affected by the 
presence of the singlet superfield. The expressions for the corresponding mass matrices are the same 
as in the MSSM, once we identify 



tan (3 



V2 
Vi 



H = Xv s 



= Xv s Ay, 



X V1V2 



(10) 



where /i is the usual Higgs mass term in the MSSM superpotential and is the corresponding term 
in the soft SUSY-breaking scalar potential. In particular, the chargino mass matrix is as in eq. (22) 
of ref. [32], and the charged-Higgs mass is 



B„ 



in 



H+ 



(11) 



sin (3 cos (3 

From eqs. (|1J), © and dSJ) it can be seen that, when A <C 1, the singlet and the singlino decouple 
from the Higgs and higgsino sectors of the MSSM, respectively. Since the experimental lower bounds 
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on the chargino masses require \i to be (at least) of the order of 100 GeV, eq. (|10p implies that, in 
this limit, v s 3> («i,«2) > so that ~ Xkv 2 is driven to rather large values - leading in turn to the 
so-called decoupling limit of the MSSM - unless k< 1 as well. 

3 One-loop corrections to the Higgs mass matrices 

In this section we describe the calculation of the one-loop corrections to the neutral Higgs boson 
masses in the NMSSM. We follow closely the approach^! of the MSSM calculation of ref. [12] . which 
is the one implemented in most public computer codes [24J that compute the mass spectrum of the 
MSSM. 

Including the one-loop corrections in the DR renormalization scheme, and using the minimization 
conditions of the scalar potential to replace the soft SUSY-breaking Higgs masses with combinations 
of the other parameters, the 3x3 mass matrices for the CP-even and CP-odd fields read 

(A4) 1 """ = (Mlf" + ± S fT,- n.,,(^) (12) 

(M^ = (Ml)"" + ^T, - U„y) (13) 

where: the tree-level mass matrices (.A/fg) tree and (.A/fp) tree are given in eqs. (|4|) and ([6j), respectively, 
and they are expressed in terms of DR-renormalized parameters; v% stands for (v±,V2,v s ); Ti is the 
finite part of the one-loop tadpole diagram for the scalar Si) n SiSj (p 2 ) and n piPi (p 2 ) are the finite parts 
of the one-loop self energies for scalars and pseudoscalars, respectively; p 2 is the external momentum 
flowing in the self energy. The explicit formulae for the scalar and pseudoscalar self energies and for 
the scalar tadpoles are collected in appendix B. We checked that, in the limit in which A — > while 
[i = Xv s is constant, our results for the 2x2 upper-left submatrix of H SiSj and for the tadpoles T\ and 
T2 coincide with the MSSM results of ref. [T2], as does the pseudoscalar self energy that we can 
obtain by rotating the 2x2 upper-left submatrix of n PiPj by an angle f3. 

The radiatively corrected squared mass of the n-th scalar or (physical) pseudoscalar can be obtained 
by solving iteratively for the n-th eigenvalue of the corresponding mass matrix evaluated at an external 
momentum p 2 equal to the mass eigenvalue itself (we remark that this procedure includes in the results 
for the masses also contributions that are formally of higher order in the perturbative expansion) . On 
the other hand, there is a well-known ambiguity in the definition of the radiatively corrected mixing 
matrices R s and R p , because the rotations that diagonalize the radiatively corrected mass matrices 
depend on the choice of external momentum in the self energies. This ambiguity reflects the fact 
that the mixing matrices themselves are not physical observables. In our analysis we will define the 

2 There are however a few differences between our conventions and those of ref. |12j : we normalize the Higgs vevs in 
such a way that (u? + «!) w (174 GeV) 2 , and our convention for the sign of the term \SH1H2 in the superpotential 
corresponds to the opposite sign of fi w.r.t. ref. [T5]. 
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radiatively corrected mixing matrix as the one that diagonalizes the mass matrix at p 2 = 0. This 
corresponds to the result obtained in the effective potential approximation. 

The tree-level mass matrices in eqs. (|12p and f)13|) depend on the combination of gauge couplings 
g 2 , on the NMSSM superpotential couplings A and k and on the three vevs v\ y vi and v s . As long 
as no experimental information on the parameters of the NMSSM Higgs sector (nor on the validity 
of the NMSSM itself) is available, A, k, v s and the ratio of vevs tan (3 can be considered directly as 
DR-renormalized inputs at some reference scale Qq. On the other hand, the DR values of v 2 = v 2 + v 2 , 
g and g' can be extracted from the experimentally known SM observables. For example, starting from 
the muon decay constant and the gauge-boson pole masses, we can make use of the relations 

v- 2 = 2V2G, ^-^^--Sy^ , (14) 

f = v~ 2 M 2 (1 + ^%^) , 9 2 = 2v~ 2 M 2 w (l+ ^^f^ ) • (15) 

In eqs. (TB|) and (fT5|) . Hy V (p 2 ) (V = denotes the finite and transverse part of the self energy 

of the vector bosons, while <5vb denotes the sum of vertex, box and wave-function-renormalization 
corrections to the muon decay amplitude. The explicit formulae for the vector-boson self energies are 
collected in appendix B. The SM contribution to <5vb was computed long ago [34], and the SUSY 
contribution can be obtained from the MSSM results given in eqs. (C.13)-(C.22) of ref. |12j . by simply 
extending the sum over the neutralinos to the five mass eigenstates of the NMSSM. 



4 Two-loop corrections in the effective potential approach 



We now discuss the computation of the two-loop corrections to the NMSSM Higgs mass matrices in the 
effective potential approach. The effective potential for the neutral Higgs sector can be decomposed 
as V e s = Vq + AV, where AV contains the radiative corrections. The 3x3 mass matrices for the 
CP-even and CP-odd fields can be decomposed as 



M 



+ ( AMI 



eff / n \ tree 

... = (Ml).. 
and the radiative corrections to the mass matrices are 



(m 



/ 1 



off 



(mI)^ + (am 



(am%) 

(AMI) 



ij 



'j 



1 

_L 5 ji 

y/2 vi 



Sij dAV 



dSi 
dAV 



dSi 



+ 



+ 



d 2 AV 



dSidSj 
d 2 AV 



dPidPj 



(16) 

(17) 
(18) 



where Vi stands for (vi,V2,v s ), and the derivatives of the correction AV are computed at the minimum 
of V e g. The comparison between eqs. (fT7]) and (fTHj) and eqs. (fT2^) and (fTHj) highlights the correspondence 
between tadpoles, self energies and derivatives of the effective potential. In the calculation of the 
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MSSM Higgs boson masses it is customary to reorganize the corrections in such a way that (A^|.) tree is 
expressed in terms of the non-zero eigenvalue of (.Mp) eff , which in the effective potential approximation 
corresponds to the physical ^4-boson mass. In the case of the NMSSM this reorganization is not as 
practical, because there are two non-zero eigenvalues of (.Mp) e . While it is possible to absorb some 
of the radiative corrections in an "effective" trilinear coupling A\, this parameter does not allow for a 
direct physical interpretation. Therefore, we refrain from this manipulation as well and leave eq. (|17j) 
as it stands. Throughout the calculation we assume that all the parameters entering both the tree-level 
and one-loop parts of the mass matrices are renormalized in the DR scheme at a renormalization scale 
that we denote by Q. 

The 0(a s ) contribution to AV from two-loop diagrams involving top, stop, gluon and gluino has 
been computed e.g. in refs. [161118] . It is the same for the MSSM and for the NMSSM, and we give it for 
completeness in appendix C. The corresponding 0{ata s ) corrections to the mass matrices in eqs. (|17|) 
and (fT8|) can in turn be computed by exploiting the Higgs-field dependence of the parameters appearing 
in AV. As detailed in ref. [18], if we neglect D-term contributions controlled by the electroweak gauge 
couplings the parameters in the top/top sector depend on the neutral Higgs fields only through two 
combinations: 



X 



\X\e iv = h t H$ , X = \X\e^ = h t (A t H$-\S*H^*^ . (19) 



The top/stop 0(a s ) contribution to AV can be expressed in terms of five field-dependent parameters, 
which can be chosen as follows. The squared top and stop masses 

1 



2 I v \2 2 

m+ = \X\ , rrir 



tl,2 2 



(m 2 Q + ml + 2 \X\ 2 ) ± yVg " ™u? + 4 \X\ 2 



(20) 



a mixing angle 9^, with < 9^ < vr/2, which diagonalizes the stop mass matrix after the stop fields 
have been redefined to make it real and symmetric 

- 2 \X\ 

sin 2 f = ' ' , (21) 

and a combination of the phases of X and X that we can choose as 

cos( y -0) = Re «' te( ^ +Im " Imm . (22) 

\X\\X\ 

A sixth parameter, the gluino mass m,g, does not depend on the Higgs background. In the following 
we will also refer to 9t, with —tt/2 < 9t < vr/2, i.e. the usual field-independent mixing angle that 
diagonalizes the stop mass matrix at the minimum of the scalar potential. 

With a lengthy but straightforward application of the chain rule for the derivatives of the effective 
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potential, the corrections to the Higgs mass matrices in eqs. (|17p and (|18p can be expressed asH 
AMl) u = l^^s^Fs + httanp ^^ F, (23) 

tl 52 

AM|) = -h\iim t s 2 e t F 2 - \h\A tl is\ et F z - h 2 f At 2 F , (24) 

tl t% 

AM 2 S ) 22 = 2 /i 2 mf Fl + 2 fc? A t m t s 2 e t F 2 + i /i 2 A| F 3 + fe? cot /? F > (25) 

51 *2 

AA"*|) 1Q = IhtXmtf^ cot(3s 2 20t F 3 - Mm /'^^ F, (26) 
' 13 2 mf — mf 

tl 52 

AJWc") = -h t \m 2 cot (3 s 2 e t F 2 - -h t XA t m t cot /3 sL. F 3 - /t t A cot /? F, (27) 

/ 23 2 mf — mf 

tl £2 

t 2\ _ 1 \2 _2 „„ + 2 « „2 Z7_ i \2 + « m t At 

i~ -n , 
ti t%' 



AM^J 33 = 2 A m * COt P F3 + COt ^ „ (m g _ m 2) F ( 2§ ) 



AM 2 P ) = h 2 tan/3 ^ At 9 F A , (29) 
/ll mf - mf 

tl 52 

A.Mp) = fr 2 2 fU, (30) 

/ 12 mf - mf 

ti t 2 

AM 2 P ) = h 2 cot/3 f At = Fa, (31) 
/22 mf -mf 

El 52 

A.Mp) = fetA 2 m< ^ 2 F A , (32) 
/ 13 mf — mf 

*1 12 

AMp) = htXcotp ™ l \ F A , (33) 
/23 mf - mf 

tl 52 

AjVip) = A 2 cot /3 . ^ At 9 . F4 , (34) 
/33 u m; -mf ) 

' v 51 52 y 

where the functions F,, F and Fa are combinations of the derivatives of AV evaluated at the minimum 
of the effective potential: 

p _ d 2 AV + d 2 AV + d 2 AV +2 d 2 AV +2 d 2 AV +2 9 2 AF 

(dm 2 ) 2 (dm 2 ) 2 (dm 2 ) 2 dm 2 dm 2 dmfdm 2 dm} dm 2 ' 

9 2 AF 9 2 AF <9 2 AV 9 2 AF 

-^2 77; 9 \o — Tci 9 \o "r 



(9m? ) 2 (9m? ) 2 dm 2 dm 2 dm 2 dm 2 

K tl' v 52 El L 52 



4c 2 6t / 9 2 AF 9 2 AF i d 2 A ^ 
' dc 2 n dm 2 dc 2 „ dm 2 dc 2 ^d 



+ fl„2 o_„2 + aJ2 fl„„2 ' (^) 



mf — mf \ dc% s omr ocz s om^ ocz s om~ 

ti 52 \ tot ti ZUt 52 ■ 



3 The differences with respect to eqs. (25)-(30) of ref. [18] have multiple origins: we do not absorb part of the corrections 
in the tree-level mass matrices; we adopt the opposite convention for the sign of /i; we take directly the derivatives of the 
renormalized effective potential as in refs. [221 \'61\ , removing the need for the counterterm- induced shifts AFi and AFi. 
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d 2 AV d 2 AV d 2 AV 2 ( dAV dAV 

+ 



(dm 2 ) 2 (dm- ) 2 dm- dm 2 m~ — m~ \ dm- dm 

v tl' x ti' tl t2 tl t2 \ tl 



16c k (2 9 2 AV „dAV\ 8c 2 et ( d 2 AV d 2 AV 



(m~ — m~ ) 2 I ^ 29t (d(% s ) 2 ^ dc 2 = ) m~ — m~ \ dc 2 ^ dm 2 z dc 2 ^ dm 2 ) ' 

ti t'2 ' \ x lot ' lot / tl 12 \ lot tl lot tl / 

p = dAV dAV 4c| 6t dAV 
dm^ dm^ m 2 — m~ dc 2 ^ ' 

tl C2 tl 12 ^ft 

F _ dAV _ dAV _ 4:c 2 9t dAV _ 2z t fi cot/? dAF 
A dm| dm| m| - m| dc^ A t ( m | _ m |) 3c^ t _0 t ' 

In eqs. (f23|) -(f39" j) above we adopted the shortcuts = cos <f> and = sin </> for a generic angle <\>. The 
parameters /i and tan/3 are defined in eq. ([TO]), and = sign(A t — fj, cot /?). 

At one loop the top and stop contributions to AV depend only on the corresponding masses. 
In units of N c /(16it 2 ), where iV c = 3 is a colour factor, the one-loop expressions for the functions 



2 

m~ 
ti 


+ m~ 

*2 


2 

ln-£ 

mj 

S2 


2 

ra- 
ti 


2 

— mj 

*2 


2 

to- 

*2 

Q 2 


-) 





appearing in eqs. (|23|) ~ (|28j) are 

22 2 
. . m~ m~ mf 

F l = ln ^4^< ^ = ln ^' ^ = 2- J 1 JM n-f, (10) 

F U = Fa = ml (log ^ " l) " < (log - 1 1 . (41) 

Inserting eqs. (|40p and (|4ip in eqs. (|23p ~ (|34p we recover the well-known results [28J for the one-loop 
top/stop corrections to the NMSSM Higgs boson masses in the effective potential approach. 

Explicit expressions for the derivatives of the contribution to AV from two-loop diagrams with top, 
stop, gluino and gluon are provided in appendix C. Rearranging the various terms, it can be shown 
that the 2x2 upper-left submatrices of AM 2 S and AM. 2 P correspond to the (D(ata s ) corrections derived 
in ref. [IB] for the MSSM in the DR renormalization scheme. On the other hand, the corrections to the 
third row and third column of the mass matrices, which are specific to the NMSSM, were not previously 
available. If the one-loop part of the corrections is expressed in terms of On-Shell (OS) parameters, 
the two-loop corrections must be supplemented with counterterm contributions that account for the 
shift from DR to OS. The required 0(a s ) shifts in the parameters mt, wi| , m f 2 > s ie t an d ^t can be 
found in appendix B of ref. [18J. 

The computation described above allows us to obtain also the two- loop 0(aba s ) corrections induced 
by the bottom/sbottom sector, which can be relevant for large values of tan/3. To this purpose, the 
substitutions t — ► b, tan/3 <-> cot/3, (AM 2 sp ) <-> [AM 2 S ^ and (AM 2 sp ) <-> {AMgp^ 
must be performed in eqs. (f23l) - (f3"9~|) . In the case of the bottom/sbottom corrections, however, passing 
from the DR to the OS scheme involves additional complications, as explained in ref. [20] . 

In the case of the MSSM, the computation of the two-loop 0(a 2 + atab + a 2 ) corrections induced 
by the Yukawa interactions of quarks, squarks, Higgs bosons and higgsinos is also available [22J. 
In contrast to the case of the 0{otta s + ai,a s ) corrections, however, this computation cannot be 
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straightforwardly extended to the NMSSM, because the Higgs and higgsino sectors are extended by 
the presence of the singlet superfield. A dedicated calculation of the 0(a^ + ata^ + a^) corrections to 
the Higgs masses in the NMSSM goes beyond the scope of this paper. 

Finally, since V e ff generates one-particle-irreducible Green's functions at vanishing external mo- 
mentum, it is clear that the effective potential approach neglects the momentum-dependent effects in 
the Higgs self energies. The complete computation of the physical masses of the CP-even and CP-odd 
Higgs bosons requires the full, momentum-dependent two-point functions (a detailed discussion of the 
correspondence between the effective potential approach and the full computation has been given in 
ref. [20J). However, in the last paper of ref. [25] it has been shown by direct calculation that, in the 
MSSM, the numerical effects of the two-loop momentum-dependent contributions to the Higgs boson 
masses are very small. There is no reason to expect that such effects would be much larger in the 
NMSSM. 

5 Numerical examples 

In this section we briefly discuss the numerical effect of the one- and two-loop corrections to the 
NMSSM Higgs masses presented in sections [3] and [H respectively. 

Among the Lagrangian parameters that enter the computation of the NMSSM Higgs masses, the 
gauge and third-family Yukawa couplings, as well as the electroweak symmetry breaking parameter 
v, can be extracted from the known values of various SM observables by taking into account the 
appropriate radiative corrections. We use the following input values for our analysis: the gauge boson 
masses M z = 91.1876 GeV and M w = 80.40 GeV; the muon decay constant = 1.16637 x 10~ 5 
GeV -2 ; the strong coupling constant a s (Mz) = 0.1189; the pole top mass Mf = 173.1 GeV; the 
running bottom mass m&(m(,) = 4.23 GeV; the tau mass m T = 1.777 GeV. Consistency with our 
computation of the one-loop radiative corrections requires that all the parameters entering the tree- 
level mass matrices be expressed in the DR renormalization scheme at a common scale Qo, which 
we take of the order of the soft SUSY-breaking scale. For consistency with the computation of the 
two-loop 0(otta s + otba s ) corrections, the top and bottom masses and Yukawa couplings entering the 
one-loop part of the corrections must also be expressed in the DR scheme. We determine the running 
electroweak gauge couplings and v directly at the scale Qo by means of eqs. (fT4"j) and (fT5|) . This 
procedure neglects the resummation of potentially large logarithms of the ratio of the weak scale to 
the SUSY-breaking scale (incidentally, we also neglect the small SUSY contributions to <5vb)> but it is 
accurate enough for the purposes of our study. The top pole mass is converted into the corresponding 
running mass, then both the top and bottom masses are evolved up to the scale Qo by means of the 
SM renormalization group (RG) equations. At that scale the SM running masses are converted into 
NMSSM running masses by the inclusion of gluino- induced threshold corrections (which are the same 
as in the MSSM). The tau mass enters only the one-loop part of the calculation and is not subject 
to QCD corrections, thus we use directly the pole mass. Finally, the strong gauge coupling a s enters 
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only the two-loop part of the calculation, therefore its precise definition amounts to a higher-order 
effect. We evolve a s from Mz to Qq by means of the SM RG equations. 

To exemplify the effect of the one- and two-loop corrections to the neutral Higgs masses in the 
NMSSM, we choose the SUSY input parameters in such a way that the scalar component of the singlet 
is relatively light and has a sizeable mixing with the lightest MSSM-like scalar. For what concerns 
the Higgs sector, we keep A as a free parameter and fix the remaining parameters as 

K = A/5 , tan/? = 2 , A x = 500 GeV , A K = -10 GeV , \i = 250 GeV , (42) 

where we take /iasa bookmark for the singlet vev v s = At/A. We adopt a common soft SUSY-breaking 
mass Ms for all of the squarks and sleptons, and fix the remaining soft SUSY-breaking parameters as 

A t = A b = A T = -1.5 M s , M 3 = 2 M s , M 2 = 2/3 M s , Mi = M s /3 . (43) 

All of the parameters in eqs. (|42p and (|43p are meant as DR running parameters at the scale Qq = Ms- 
Figs. [1] and [2] exemplify the effect of the one-loop corrections to the NMSSM scalar masses. In 
fig. [1] we plot the squared rotation matrix elements (-R13) 2 and (-Rfs) 2 ' w hich measure the strength of 
the singlet component in the two lightest scalars hi and h 2 , as a function of A for Ms = 300 GeV. 
In fig. [2] we plot the masses of the two lightest scalars for the same choices of inputs. In both plots, 
the dotted lines correspond to the tree-level results; the dashed lines include the one-loop 0{at) and 
G(atb) corrections computed in the effective potential approach; finally, the solid lines correspond to 
the results of the full one-loop calculation. For the full one-loop calculation of the rotation matrix the 
external momentum in the scalar self energies is set to zero. It can be seen in fig. [I] that, at small A, 
the lightest scalar hi is dominantly MSSM-like while h 2 is dominantly singlet. When A increases the 
mixing between singlet and lightest MSSM-like Higgs increases as well. Meanwhile, the heaviest scalar 
A13 has a mass of the order of 600 GeV and its singlet component is always small. It is interesting to 
note that - at least in this point of the parameter space - the value of A for which the two lightest mass 
eigenstates cross over (i.e., hi becomes dominantly singlet) depends quite strongly on the accuracy of 
the calculation. In particular, when only the quark/squark contributions to the radiative corrections 
are included the crossover occurs for much lower values of A than in the full one-loop calculation. 
Fig. [2] shows the effect of the radiative corrections to the two lightest scalar masses. The rise with A 
in the tree-level masses is due to the well-known NMSSM contribution to the Higgs quartic coupling 
proportional to A 2 sin 2 2(3. The comparison between the dotted and dashed lines shows that the G(ctt) 
corrections induced by top and stop loops have a particularly large effect on for small values of A, 
when hi is light and mostly MSSM-like. The 0(ab) corrections induced by bottom and sbottom loops 
are also included in the dashed lines, but they are negligible due to the small value of tan (3. When A 
increases the 0{att) corrections are shared between m/ ll and m^ 2 , and become less relevant due to the 
increase in the tree-level masses. However, even for the heavier scalar h 2 these corrections can still 
amount to several GeV at large A. Finally, the comparison between the solid and dashed lines in fig. [2] 
shows that the remaining one-loop corrections - which constitute one of the original contributions of 
this paper - are also relevant, and can account for shifts of 5-10 GeV in both masses. 
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Figure 1: The squared rotation matrix element (Rf^) 2 , measuring the singlet component in the scalars 
h\ and h,2, as a function of A, for Ms = 300 GeV. The values of the other input parameters and the 
meaning of the different curves are described in the text. 
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Figure 2: The masses of the two lightest scalars /ii and ^2 as a function of A, for Ms = 300 GeV. The 
values of the other input parameters and the meaning of the different curves are described in the text. 
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Figure 3: The squared rotation matrix element (R? 3 ) 2 , measuring the singlet component in the scalars 
hi and h%, as a function of Ms, for A = 0.5. The values of the other input parameters and the meaning 
of the different curves are described in the text. 
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Figure 4: The masses of the two lightest scalars hi and hi as a function of Ms, for A = 0.5. The 
values of the other input parameters and the meaning of the different curves are described in the text. 
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Figs. [3] and 0] exemplify the effect of the two-loop corrections to the NMSSM scalar masses. In 
fig. [3] we plot (-R13) 2 and (-R23) 2 as a function of M$ for A = 0.5. In fig. [5] we plot the masses of the 
two lightest scalars for the same choices of inputs. In both plots, the dotted lines corresponds to the 
full one-loop results (again, the rotation matrix is computed at zero external momentum); the dashed 
lines include the two-loop leading-logarithmic 0{ata s ) contribution to the (2,2) entry of the scalar 
mass matrix as implemented in NMHDECAY |30j . i.e. 



Finally, the solid lines correspond to the results of our two- loop 0{pna s + ai,a s ) calculation. It can be 
seen in fig. [3] that for small Ms the lightest scalar h\ is mostly MSSM-like while h,2 is mostly singlet. 
When Ms increases, the radiative corrections increase the mixing between singlet and MSSM-like 
Higgs. Fig. U] shows that the two-loop corrections to the lightest scalar mass are positive and relatively 
small. This is a typical feature of the DR computation, in contrast to the OS computation in which 
the two-loop corrections are negative and much larger (for a discussion of this issue in the MSSM see 
ref. [35J). It is interesting to note that, in this scenario, the leading-logarithmic term accounts only for 
a fraction (30% to 60%, increasing with Ms) of the total 0(ata s ) contribution to the (2,2) entry of 
the scalar mass matrix. Indeed, the leading-logarithmic approximation of eq. (|44|) neglects potentially 
large contributions controlled by powers of the ratio At /Ms, as well as the possibility of mass splittings 
among stops and gluino. The effect of the 0(ata s ) corrections to the entries of the scalar mass matrix 
other than (2,2) is also non-negligible. The comparison between the dashed and solid curves for hi in 
figs. [3] and H] shows that, in this point of the parameter space, the non-leading-logarithmic contributions 
contained in our two-loop calculation induce a shift of 1-2 GeV in m/ ll , and have a sizeable effect on 
the mixing matrix as well (on the other hand, the near overlap of the dashed and solid curves for /12 in 
fig. 0]is the result of an accidental cancellation). One of the attractive features of the NMSSM is the 
viability of scenarios in which Ms is not much above the weak scale. It is clear from figs. [3] and H] that, 
in those scenarios, the leading-logarithmic approximation is not satisfactory, and a reliable evaluation 
of the two-loop corrections requires at least the complete 0(ata s + abtt s ) calculation. 

To conclude this section, we show in fig. [5] the effect of the radiative corrections to the mass of 
the lightest physical pseudoscalar A\. For the choice of parameters considered in this example A\ 
is almost entirely singlet, therefore its mass is hardly affected by the one- and two-loop corrections 
involving quark/squark loops. The pseudoscalar A2, on the other hand, is almost entirely MSSM-like, 
but its tree-level mass is of the order of 600 GeV, thus it is also not much affected by the radiative 
corrections. However, the Higgs self-interactions and the Higgs-higgsino interactions controlled by 
the superpotential couplings A and n do induce non-negligible corrections to the lightest pseudoscalar 
mass. The dashed and solid lines in fig. [5] correspond to the tree-level and one-loop determinations of 
mA 1 , respectively, as a function of A. The input parameters are chosen as in eqs. (|4"2"j) and (@3|), but 
we show two sets of curves corresponding to k = A/5 and k = A/3. From the comparison between the 
dashed and solid curves it can be seen that the one-loop corrections to tua 1 can amount to several GeV 




(44) 
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Figure 5: The masses of the lightest physical pseudoscalar A\ as a function of A, for Ms = 300 GeV 
and k set equal to either A/5 or A/3. The values of the other input parameters and the meaning of 
the different curves are described in the text. 

when A and k take on relatively large values. The two-loop corrections computed in this paper include 
only the quark/squark contributions, therefore the corresponding curves would essentially overlap with 
the one-loop curves. 

6 Conclusions 

The NMSSM is an attractive extension of the MSSM: it provides an elegant solution to the \i problem, 
it reduces the need for heavy superpartners to lift the Higgs mass through radiative corrections and it 
has an interesting collider phenomenology. However, the accuracy of the theoretical predictions for the 
NMSSM Higgs masses has until now been stuck to the level that for the MSSM had been achieved in the 
mid-1990s. In this paper we took a few steps towards bridging the accuracy gap between the NMSSM 
and MSSM calculations. In particular, we provided a full one-loop calculation of the self energies and 
tadpoles of the neutral Higgs bosons of the NMSSM, and we computed the two-loop 0{a,%a s + otba s ) 
corrections to the neutral Higgs boson masses in the effective potential approximation. We showed 
that both classes of corrections can induce shifts of a few GeV in the light scalar and pseudoscalar 
masses, and they can also sizeably affect the mixing between singlet and MSSM-like Higgs scalars. 
Taking these corrections into account in phenomenological analyses of the NMSSM Higgs sector will 
be crucial for a meaningful comparison between the MSSM and NMSSM predictions. 
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Appendix A: definitions of the couplings 

In this appendix we provide explicit formulae for the couplings that enter the calculation of the one- 
loop corrections to the Higgs boson masses in the NMSSM. 

Higgs sfermion couplings: the terms in the NMSSM Lagrangian relevant to the calculation of 
the sfermion contributions to the Higgs self energies can be written as 

c D -Y. X s lSj F k F e SiSjF^Fe -^2\ iPj p k F e p i p J F k F e -Y. X Sl F k F e S i F k F l - Y. iX Pl F k F e P i F k F i , 

ijkl ijki iki ikl 

(Al) 

where F{ = represent the sfermions in the basis of interaction eigenstates. The quartic 

couplings in eq. (|Aip are symmetric with respect to the exchange of i and j as well as with respect 
to the exchange of k and 1. The trilinear couplings of the scalars are symmetric with respect to the 
exchange of k and i, whereas the trilinear couplings of the pseudoscalars are antisymmetric. The 
couplings to up-type and down-type squarks (the generalization to the sleptons is straightforward) are 

X - X l! \ - X - 

A sisiC/i(7i ~~ A pipiUiUi ~ 2 Sul ' S1S1U2U2 ~ A pipiU 2 U 2 ~ 2 9ur ' 

X s 2 s 2 U 1 U 1 ~ \ 2 p 2 UiUi ~ ~~^9u L + , A s 2 s 2 U 2 U 2 ~ X p 2 p 2 U 2 U 2 ~ 2 9ur + ~2~ ' 

X sis z U 1 U 2 = ~ \iP S UiU 2 = T~ ' ( A2 ) 



A s 1 s 1 D 1 D 1 ~ A p 1 p 1 D 1 D 1 ~ ~~2~ 9d L + -n- j A s 1 s 1 D 2 D 2 ~ A pipiD 2 D 2 ~ ~^ 9d R + -~~ J 



X X - X X -t- 

A s 2 s 2 b 1 bx — A p 2 p 2 b x bx — 2 9dL ' A s 2S2 b 2 b 2 ~ A P2P2 b 2 D 2 ~ 2 9dR ' 

A s 2 s 3 b 1 D2 = ~ \2p3D1D2 = 4~ ' ( A3 ) 

X si&iUi = ^9 2 9u L v u X . Sl u 2 U 2 = ^9 2 9u R vi, 
X s 2 u 1 u 1 = -^9 2 9u L v 2 + V2hlv 2 , x S2 u 2 u 2 = -^9 2 9u R v 2 + V2h 2 u v 2 , 
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\ - \ - hu Vs A \ - \ - hu Au \ - \ - _huV\\_ 

A s 1 u 1 u 2 - A Pl u 2 Ui - k ' A s 2 u l u 2 ~ a p 2 vJj 2 - k > A s 3 ^iC/2 ~~ \zV 2 Vi - ^ ' 

(A4) 

X s 1 d 1 d 1 = V2f9d L vi + V2h%vi, ^ Sl b 2 b 2 = V2f 9d R vi + V2h 2 d v 1} 

X a 2 D 1 Di = -^9 2 9d L v 2 , K 2 b 2 b 2 = -^9 2 9d R v 2 , 

\ - \ hd Ad \ - \ hd VsX \ - \ - hd V2A 

A s 1 D 1 D 2 - A p 1 D 1 D 2 - ' A s 2 D 1 D 2 - A p 2 £>2.Di ~~ /2 ' A s 3 Dl£>2 ~ A p 3 D 2 Di ~ K ' 

(A5) 

where gt = J3 — et sin 2 0\y. Here I3 is the weak isospin and e/ is the electric charge of the chiral 
superfield that contains the sfermion (e.g., e UR = —2/3). The couplings that cannot be obtained 
by swapping the last two indices and - for the quartic couplings - the first two indices of those in 
eqs. (|A2p ^ (|A5p are all vanishing. 

In the absence of CP-violating phases in the sfermion mass matrix, the sfermion mass eigestates 
fi = (/11/2) are related to the interaction eigenstates Fi by an orthogonal rotation: 

\ —sin Of cos Of J 

where Of is the mixing angle of the sfermions /. The quartic and trilinear couplings between Higgs 
fields and sfermion mass eigenstates are related to the corresponding couplings between Higgs fields 
and sfermion interaction eigenstates as follows 

A 4>i4>jhh = R ka R ®> A <t>i<t>jF a F b ' \ifkft = R ka R £b X <f> t F a F b ' ( A7 ) 

where <j>{ represents either Sj or pi, and summation over repeated indices is understood. 

Higgs self-couplings: the terms in the NMSSM Lagrangian relevant to the calculation of the 
neutral-Higgs contributions to the Higgs self energies can be written as 

C D — ^ A SiSjSfcSf SiSjSkSi — ^PiPjp k pi PiPjPkPi — X] As i s oPkPi SiSjPkPe 

ijki ijki ijki 

~ ^2 x s iSk s e SiS k Se - ^2 x s iPkPt SiP k Pe ■ (A8) 

iU ikl 

The quartic and trilinear neutral-Higgs self couplings entering eq. (|A8|) are symmetric with respect to 
the permutation of any two indices corresponding to fields of the same parity. They read 

A — A A -A - t- A — A - — 

AsisiSISJ — ^S 2 S 2 S 2 S 2 — ^pipipipi — ^p 2 p 2 p 2 p 2 — 1 A S3S3S3S3 — A P3P3P3P3 ~ ^ > 



A A - — (2\ 2 -n 2 \ A — A - — 

A s 1 s 1 s 2 s 2 — A pypip 2 p 2 — i/ J > / ^s 1 s 2 S3S3 — A pip2P3P3 ~~ 24 ' 
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A 2 

<^SlSlS3S3 = A S2S2S3S3 = ^P1P1P3P3 = A P2P2P3P3 = ^ ' 

X X - Q— X - X -I /Sx 2 -n 2 

^SiSipipi — A S 2 S 2 p2P2 — a ' A SlS!P2P2 — A S 2 S2PlPl o I y 



^SlSlp 3 p 3 — ^S2«2P3P3 ~~ -^S 3 S 3 pipi — A. 



S3S3P2P2 



-^SiS 2 P3P3 ~~ ^«3S3PlP2 ~~ -^SiS 3 p2P3 ~ A S2S3PlP3 ~ A ! -^S3«3P3P3 ~~ ! (^-9) 



f «1 



AsjsiSl — -^Sipipi — n /o ' ^S 2 S2S2 — A S2P2P2 



2^2 

A SVP2P2 ~ A SiS2S2 — 2yJ~2 V / ' s 2PlPl ~~ "^SaSlSl — 2^/2 \ 

A 2 « 



1 

8 


(2A 2 - 


A 2 




4 


5 


-^S3S3P3P3 








5 2 ^2 




2\/2 




W2 



Ak k 2 

"J" ' -^S 3 S3P3P3 = ' 



•^S3pipi — ^S3P2P2 ~~ ^ -^S 3 SiSi — 3 A S j 



3«2S2 



x/2 ' 

As 3 S3S3 = g^/^ + V 2 K t> s , A S3 p 3 p 3 = — ^" + ViK 



'Vs ^ ^=(AVl-KU 2 ) , A,,,, ,,,, = ^= (A Co- ) 



x A A . . Aki> 2 Akui 

-^«1P3P3 ~~ "^/^ {AV1+KV2) , A S2 p 3 p3 — — J= {AV2 + KV1) , A S3 p ip3 — — y= , A S3 p 2 p 3 — — ^= , 

x A A A A k v s _ \A\ Xkv s _ \A\ | Xkv s 

S1S2S3 ~ 6V2 3V2 ' S1P2P3 ~ S2P1P3 " 2V2 V2 ' S3P1P2 " 2^ + V2 ' 

(A10) 

All of the couplings that cannot be obtained by permuting the indices of the couplings in eqs. (jA9j) 
and (|A10p vanish. Rotating two scalar or pseudoscalar interaction eigenstates into mass eigenstates as 
in eqs. ([5]) and (|7j), and exploiting the permutation symmetry of the original couplings, the couplings 
that enter the calculation of the scalar self energies can be expressed as 

A SiSjh k h e = 6 Rgfr \ SiS:jSaSb , A SiSj(L k ai = Rka P lb A SiSjp a Pb > 

A Sih k h t = 3 R% a Rji b \ SiSa s b , A Si a k a e = ^ka ^ib A SiPaPb ' (All) 

Similarly, the couplings that enter the calculation of the pseudoscalar self energies are 

A PiPja k a e = 6 Rka-^lb A PiPiPaPb ' A PiPjh k h e = Rka P £b A s a s b PiPj > 

A Pia k h e = 2 Rf a R% b \ Sa p b pi • (A12) 

The terms in the NMSSM Lagrangian relevant to the calculation of the charged-Higgs contributions 
to the Higgs self energies can be written as 

^ D ~22K i s j h+hJ S i S jKh ~22\iPjh+hJ P i P jK h £ 
ijkl ijk£ 

- E Kh+h- so%hj - £ i x Pihth - ^ , (A13) 

iM ikl 
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where we express the charged Higgs fields directly in the basis of mass eigenstates, so stands for 
(G^ , H^). In our conventions the relation between the charged Higgs mass eigenstates and interaction 
eigenstates reads @ 

V H ± ) \ Hf J \ sin (5 cos (5 ) 

The quartic couplings in eq. (|A13|) are symmetric with respect to the exchange of i and j as well as with 
respect to the exchange of k and £. The trilinear couplings of the scalars are symmetric with respect 
to the exchange of k and £, whereas the trilinear couplings of the pseudoscalars are antisymmetric. 
Explicit expressions for the couplings are 

X s lSl h+h- = X s 2 s 2 h+h- = X PlP1 h+h- = X P 2P2h+h~ = g (9 2 + g' 2 cos2(5) , 

Ksih+h~ = X sosoh+h7 = Kim h+hZ = Kovih+h.7 = o ~ 9 2 COS 2(5) , 



S2S2h^h 1 ' p\p\h 2 h 2 " P2P2h^h 1 

X s 1S2 h+h- = X pi P 2h+h- = ~ X s 1S 2h+h- = ~ X PiP 2h+h- = „(2A 2 -3 2 )sin2/? 



g' 2 

A u+u— — A 1.+7,— — — A u+u— — — A u+u— — — sin2/3 , 

sisih~h 2 p\p\hih 2 S2S2nJh 2 P2P2h{h 2 g • ' 

Ks 2 h+h- = - X p, P 2hih- = \( 2x2 -9 2 ) cos 2/? , X S3S3h + h - = - x P3P3h + h - = cos 2/3, 

= ^3^+^ = ^ (A sin 2/3) , A S3S3ft + fe - = A^^- = ^ (A + k sin 2/3) , (A15) 

1 / ' 2 i „/2„„„r,<A , ,. /o\2 „2\ 



X s\h 2 h 2 



Khih- = -^m (vi(g' + g' 2 cos2P) + v 2 (2X 2 - g') sin 2/3 
^ (^i(5 2 -5 /2 cos2/3) - W2 (2A 2 - 9 2 ) sin 2/3 
= ^ (-^i5' 2 sin2/3 + V2 (2A 2 - 5 2 )cos2/3) , 
\^u+u- = (v 2 (g 2 - 5 ' 2 cos2/3) + wi(2A 2 - 9 2 ) sin2/3 



^= (v2 (g 2 + g' 2 cos 2/3) - Wl (2A 2 - 9 2 ) sin 2/3) , 
= ^= (W 2 sin 2/3 + Wl (2A 2 - 5 2 ) cos 2/3) , 
x s 3h + h - = ~j= (2Av s - {A x + 2kv s ) sin 2^) , A Ssh +^ = -^= (2A v s + (^4a + 2k v s ) sin 2(5) , 

A s 3 h+h- = " (^a + 2k «,) cos 2(5 , 



4 This differs from the conventions of Ref. [12] by one field redefinition. 
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A . + .- = -^=(2\ z -g 2 ), A ,+,- = -^=(2\ 2 -g 2 ), A .+._ = -^={A x -2kv s ) . (A16) 

All of the couplings that cannot be obtained by permuting the indices of the couplings in eqs. (|A15P 
and (|A16P vanish. 

Higgs neutralino couplings: In the formalism of two component spinors, the terms in the NMSSM 
Lagrangian relevant to the calculation of the neutralino contributions to the Higgs self energies can 
be written as 

£ 3 - S i ~ i \ii>i*l>l Pi V'fcV'/ + hx - ' ( Al7 ) 

ikl ikl 

where V? = ( — ib, —iw°, h®, h®, s) are the neutralino interaction eigenstates. The couplings in 
eq. (|A17|) are symmetric with respect to the exchange of the neutralino indices k and £, and read 

K^^l = ~ ^82^1 = - X pii>^° 3 = \>2^^ = f ' ^S 3 ip0^0 = A^^O^O = ~J= , 

All of the couplings that cannot be obtained by permuting the neutralino indices of the couplings in 
the equation above vanish. The couplings of the Higgs bosons to the neutralino mass eigenstates Xi 
are related to the corresponding couplings to the neutralino interaction eigenstates as follows 

X <t>iX k X° e = N ka N Ib ) ( A19 ) 

where fa represents either Sj or pi, N is the rotation matrix defined in eq. ([9]) and summation over 
repeated indices is understood. In the absence of CP-violating phases in the neutralino mass matrix, 
the couplings can be taken as real if we allow for negative neutralino masses. 

Higgs chargino couplings: The terms in the NMSSM Lagrangian relevant to the calculation of 
the chargino contributions to the Higgs self energies can be written as 

c => -E A Si ^^- witi - j EW Pi ^J + h - c - > ( A2 °) 

ike ike 

where ipf = (—iw + , J12) and ip~ = (— iw~, hi) are the positive and negative chargino interaction 
eigenstates, respectively, in the formalism of two-component spinors. The only non-zero couplings in 
eq. (fA~20j) are 

\ _ \ _ x _ _ \ _9_ 

X s s ^+^- = \ 3 i,+4>- = ^ • ( A21 ) 
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The mass matrix M x for the charginos can be diagonalized by a biunitary transformation 

diag(m Y+ ) = U* M X V^ , (A22) 

where the unitary matrices U and V rotate the negative and positive chargino interaction eigenstates, 
respectively, into the corresponding mass eigenstates 

XT = Uijiij , xt = Vijrf . (A23) 

The couplings of the Higgs bosons to the chargino mass eigenstates are related to the corresponding 
couplings to the chargino interaction eigenstates as follows 

Wx7=^^W- (A24) 
where <f>i represents either Sj or pi and summation over repeated indices is understood. In the absence 
of CP-violating phases in the chargino mass matrix, the couplings can be taken as real if we allow for 
negative chargino masses. 

Z— neutralino and Z chargino couplings: Since the singlino is neutral with respect to the MSSM 
gauge sector, the couplings of the Z boson to the neutralinos (and, of course, to the charginos) in the 
interaction basis are the same as in the MSSM. In the formalism of two-component spinors they read 

c d -^(a^ ? ^^° + h.c.) Ev*; z ^>^ zZ x z^jZ^t^j ■ 

ij ij ij 

(A25) 

The only non-zero couplings are 

- ^ g 

x z^° = ~ X z^\ = X zi,+ip+ = X z^~ = -j^-g-i X zf+f+ = X z^,-^ = ^ cos2^. 

(A26) 

In terms of the neutralino and chargino mass eigenstates, the couplings read 

X z x »x<] = N ikN* e \ z ^° , X z x + X + = v ikV* e X Zi ,+^+ , X Zx : x - = UikUjt\ z ^ T • (A27) 

W— neutralino chargino couplings: The couplings of the W boson with charginos and neutralinos 
are the same as in the MSSM, but we give them for completeness. In the formalism of two-component 
spinors, they read 

£ -Ev^^^^° -EV 7 ^^^° + h - c - ( A28 ) 

ij ij 

The only non-zero couplings are 

X W4>°i>+ = X Wijj°^~ = ~ 9 , X Wii>°ip+ = ~ X w^- = ^ , (A29) 
In terms of the neutralino and chargino mass eigenstates, the couplings read 

Vx l >x+ = N *k v ji x w^+ > X w x °x~ = N ^ u Ie x z^~ • ( A3 °) 
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Appendix B: one- loop self energies and tadpoles 



In this appendix we list the explicit formulae for the one-loop self energies and tadpole diagrams that 
are necessary to the calculation of the neutral Higgs boson masses. The calculation is performed in 
the 't Hooft-Feynman gauge, in which the Goldstone bosons and the ghosts have the same masses as 
the corresponding gauge bosons. 



Scalar self energies The contributions to the scalar self energies from matter-fermion loops and 
from loops involving gauge bosons or ghosts are essentially the same as in the MSSM, and read [12] 

167r 2 n£J> 2 ) = 3h 2 b Up 2 -4m 2 )B (m b ,m b )-2A (m b ) 



+ 



+ 



h 2 T [(p 2 -Am 2 T )B (m T ,m T )-2A (m T ) 



E {R^f F(m h± ,M w ) + 9 - ]T (R^f F{m ail M z ) 



i=l 



i=l 



• c} [g 2 B (M W , M w ) + g 2 M 2 Z B (M Z , M z ) 



+ 2g 2 A (M w ) + 2g 2 A (M z ) 



(Bl) 



16 vr 2 Uf' 2 Y 2 (p 2 ) = 3h 2 [{p 2 -Am 2 ) B (m t ,m t ) -2A Q {m t 



+ T E F(m h ±,M w ) + V E F(m ai ,M z ) 



i=l 



i=l 



[g 2 M^y B (M W , M w ) + g 2 M 2 Z B (M Z , M z ) 



2g 2 A (M w ) + 2g 2 A (M z ) 



(B2) 



16 vr 2 !!^ 2 



i=l 



i=l 



7 

+ - s p cp [g 2 M 2 , B (M W ,M W ) + f M 2 B (M Z ,M Z ) 



(B3) 



where we neglected the Yukawa couplings of the first two generations. The loop functions Aq^ui), 
Bo(rrii,m2) and -F(mi, 7712) are defined in appendix B of ref. [12] . and they depend also on the external 
momentum p 2 and on the renormalization scale Q. Here and in the following we use for brevity the 
notation = cos0, = sin0, for a generic angle 4>. Also, hf stands for (G^, H^), a, stands for 
(G°, A\, A2), and the rotation matrices R° and R p are defined in eqs. (|A14p and (|7|), respectively. The 
contributions to the self energies involving the field £3 are zero because the singlet does not couple to 
the gauge sector nor to the matter fermions. 
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The sfermion contributions to the scalar self energies read 



16vr 2 n/ Sj V) = E E 2N cK lS] f k f k M^f k ) + E E N iKjJ t KjJ k Mrn f ~ k ,m h ), (B4) 

/ k=l f k,l=l 

where the first sum in each term runs over all the sfermion species and the second over the sfermion 
mass eigenstates. N[ is the number of colours for the sfermions /, while nj is 1 for the sneutrinos 
and 2 for all of the other sfermions. The quartic and trilinear Higgs-sfermion couplings are obtained 
by combining eqs. (lA2[) - (|A7j) . 

The Higgs contributions to the scalar self energies read 



16vr 2 n^(p 2 ) = E 2 Ks J h k h k A {m hk )+ E 2 Kh k h l K J h k h l B {m hk ,m he 



k=l 



kl=l 



+ E 2 ^s lSj a k a k A (m ak ) + E 2)[ 's l a k a e Xs :j a k a e Bo(ma k ,m ae ) 
k=l k,l=l 
2 2 

E 2X s iSj hih-Mm h ±) + E K.h+hj K.h+h; B o(m h ±,m h ±) , (B5) 



n=l 



m,n=l 



where the neutral couplings \ SiSi <t> k <t> t and \ Si <f> k <f>i, where <j) represents either h or a, are obtained by 
combining eqs. (|A9p - (|All|) . while the charged couplings X s s h+h ~ and X s . h + h - are given in eqs. (|A15P 
and (|A16p . respectively. 

Finally, the chargino and neutralino contributions to the scalar self energies read 



16vr 2 n* V 



5 

4 E 

kl=l 



+ 2 E 

kl=l 



Re ^x° fc x? X s>x° k x°J G ( m x° k > m x° e ) 



2 m^o m v o Re(A„ ^o^o A„ v o v o) Bnim^o , m^o) 

X k Xi v b iX k Xi b ]X k Xi' " v X k Xi' 



Re (K v+v- X s ■Y+y-) G ( m Y ± ' m y ± ) 

s iX k Xe h iX k Xi x k Xi 



2m xt m xf R <Kxtx7 X s jX UZ )B ^ m x^ m xt) 



(B6) 



where the loop function G(mi, m-i) is defined in appendix B of ref. [12], the Higgs-neutralino couplings 
are obtained by combining eqs. (1A18I) and (|A19j) . and the Higgs-chargino couplings are obtained by 
combining eqs. (]A21j) and (IA24I) . 

Pseudoscalar self energies The contributions to the pseudoscalar self energies from matter- 
fermion loops and from loops involving gauge bosons or ghosts read 



i6 7r 2 n££> 2 



3 h 2 b 



p 2 B (m b ,m b ) - 2A (m b ) + h 2 T p 2 B (m T ,m T ) - 2 A Q (m 
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+ y E («fi) 2 ^K±,M W ) + y E (^fi) 2 ^.M,) 

j=l i=l 

a 2 

+ \<%M%r B (M W , M w ) + 2g 2 A (M w ) + 2g 2 A (M z ) , (B7) 



16^ 2 Il£J 2 (p 2 ) = 3h 2 [p 2 B (m t ,m t )-2A (m t ) 



+ y E ( fi S) 2 F (m ht ,M w ) + ^ E (^) 2 ^,M Z ) 

j=l i=l 



O 2 

+ |-4M^ J B (iV%,M W /) + 2 5 2 Ao(M w ) + 25 2 Ao(M z ) , (B8) 

a 2 2 a 2 3 

167r 2 n££> 2 ) = y -Y,Z° x I$ 2 F{m h ±,M w )- y -Y,^2nrn hi ,M z ) 

i=l i=l 

a 2 

- y -c p spMlrB (M w ,M w ) , (B9) 

the rotation matrices R c and i? s are defined in eqs. (jA14j) and ©, respectively. The contributions to 
the self energies involving the field P3 are zero because the singlet does not couple to the gauge sector 
nor to the matter fermions. 

The sfermion contributions to the pseudoscalar self energies read 

nj ^ nj 

167r 2 n(„ (p 2 ) = VV 2JV/A f f A (m f ) - V V JV/ A f f A f f B (m f ,m f ) , 
/ fc=l / fc,£=l 

(BIO) 

where the quartic and trilinear Higgs-sfermion couplings are obtained by combining eqs. (|A2|) - (|A7p . 
The Higgs contributions to the pseudoscalar self energies read 

3 3 



16vr 2 n^ pj (p 2 ) = E 2X PiPjh k h h A (m hk ) + E 2 K P] a k a k M m a k ) 

k=l k=l 
3 

+ E \ l a k h t \ J a k h l BQ{m ak ,m ht ) 
k,e=i 

2 2 

+ E 2X p iPj h+h-Mm h ±) ~ E \ ih +h7 \ jh +h- B o(™ h ±,rn h ±) , (Bll) 



"k'H 

n=l m,n=l 



where the neutral couplings \ PiPj h k h e , X PiPjakat and Xp iak h e are obtained by combining eqs. (|A9|) and 
(lAlOh with eq. (IA12I) . while the charged couplings A„ _ h + h - and A„ ,+, - are given in eqs. (IA15I) and 

PiPj' l k Pi k H 

(|A16|) . respectively. 

Finally, the chargino and neutralino contributions to the pseudoscalar self energies read 



16vr 2 nx pj (p2) = 4 E 
ki=i 
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2 E 



kl=l 



+ 2m v o m v o Re(A„ v o v o A„ v o v o) -Bo (m v o , m v o ) 



+ 2m ± m + Re(A + - A + - ) B (m ± , m +) 



(B12) 



where the Higgs-neutralino couplings are obtained by combining eqs. (IA18D and (IA19D . and the Higgs- 
chargino couplings are obtained by combining eqs. (|A21j) and (IA24D . 

Scalar tadpoles The contributions to the scalar tadpoles from matter-fermion loops and from loops 
involving gauge bosons or ghosts are the same as in the MSSM, and read [12] 

16^ 2 T/' y = -6V2h b m b A {m b )-2V2h T m T Ao(m T ) + ^(g 2 A (M w ) + g 2 A (M z )) , 

(B13) 

16vr 2 r/' y = -QV2h t m t A (m t ) + ^[g 2 A (M w )+g 2 A (M z )) . (B14) 

The contributions to T3 are zero because the singlet does not couple to the gauge sector nor to the 
matter fermions. 

The contributions to the tadpoles from loops involving sfermions or Higgs bosons read 



Wn 2 T/>* = E E N ^Kf k f k Mm k ) + E E M^ k ) , (B15) 



k=l d> k=i 



where the first sum in the second term runs over all the Higgs bosons (i.e., (f) = h,a,h^), and 
is 2 for and 3 for h and a. The trilinear Higgs-sfermion couplings are obtained by combining 
eqs. (|A4|) - (|A7P ; the trilinear neutral-Higgs self-couplings are obtained by combining eqs. (|A10|) and 
(|A11|) ; the trilinear couplings with the charged Higgs bosons are given in eq. (|A16j) . 

Finally, the contributions to the tadpoles from loops involving neutralinos or charginos read 

5 2 
16 vr 2 Tf = - 4 £ A SjX o x o m^o io(™ x o) - 4 E K lX + Xk m x± M™^) , ( B1 6) 

k=l k=l 

where the Higgs-neutralino couplings are obtained by combining eqs. (|A18|) and ()A19p . and the Higgs- 
chargino couplings are obtained by combining eqs. (|A21|) and (jA24j). 

Vector boson self energies Since the singlet superfield is neutral under the MSSM gauge group, 
the vector boson self energies in the NMSSM differ from the corresponding MSSM quantities only 
through the effect of the mixing of the singlet and singlino with the MSSM Higgs bosons and higgsinos, 
respectively. The transverse parts of the W and Z self energies appearing in eqs. (|14H and (|15j) can 
be obtained by inserting in eqs. (D.4) and (D.9) of ref. |12j appropriate combinations of the Higgs 
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scalar and pseudoscalar mixing matrices, as well as the chargino and neutralino couplings denned in 
appendix A. 

The transverse part of the Z-boson self energy reads 



16 tt 2 U T zz (p 2 ) = 2 g 2 M 2 J2 (Rfi + Ri2 spf B (m hi , M z ) 

i=l 

3 2 

2? E (r^R^-R^R^) 2 B 22 (m hi ,m aj ) - 2g 2 c 2 9w ^B 22 i 



- Afci [2p 2 + M^-M 2 z ^\ B (M W ,M W ) - 16 f A B 22 (M W , M w ) 



r2 S w 



2 A 



\ 9w / 

2 9 2 E N c +9} R ) H(m f ,m f ) -4g fL g fR m 2 f B (m f ,m f ) 

f 

n f 2 

8 9 2 E E N c (gf^iRji-gf^R^) B 22 (m fi ,m f .) 

f ij=l v J 



4 E I^Zx°x°| 2 H(m x o,m x o) — 2m x o m x q Bo(m x o,rn x o) 
2 



+ E 



H(m ±,m ±] 



4 m + m + Re A* + + A 



B (m ±,m ±; 

At A j 



(B17) 



where the loop functions B 22 (mi,m 2 ) and H(mi,m 2 ) are defined in appendix B of ref. [12], is 
defined after eq. (1A5I) . and the couplings A^ v o v o, A 7 + + and X 7 - - are obtained by combining 

i 3 X-j Xj Xi Xj 

eqs. (TX261) and (|X27|) . 

Finally, the transverse part of the VF-boson self energy reads 



16 TT 2 l(^ W {p 2y 



= g 2 M$ v E (Rfi cp + Ri2 spf B (m ht ,M w ) 
i=l 

3 2 2 

" EE {RliR^i-R^Rh) B 22 (m hv m h ±) 

i=l j=\ 3 

~ 5 2 EE (^fl + ^2 2 ^22 (m ai , m h ± ) 
i=i j=l 3 

g 2 [(4p 2 + M 2 , + M 2 ) c 2 ew - M 2 sj w ] B Q (M Z , M w ) 
- 8g 2 c 2 9w B 2 2(M z ,M w ) - g 2 s 2 w [s B 22 {M W , 0) + Ap 2 B (M W , 0) 

3^2 ^2 

+ -5- E H ( m u,rn d ) + — E-^( m e>°) 



u/d 



e/v 
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- 6s 2 E E (^i^) 2 ^K.^) - V E E (^l) 2 B 22( 

u/d *iJ =1 i=l 

5 2 r ( \ 
+ E E [ V x ?xtl 2 + lVx«x-l 2 J #K°' m xf ) 



me., mo) 



i=l j=l 



+ 4m y o m y ± Re ( A^o x + Vx°x" ) B o{m x o,m y ± 



(B18) 



where the sums in the fermion contributions and the first sums in the sfermion contributions run over 
the three families of (s)quarks and (s)leptons, and the couplings ^w\°x + ana " ^VFx°x~ are obtained by 
combining eqs. (|A29|) and (|A30p . 



Appendix C: derivatives of the two-loop effective potential 

We provide in this appendix the explicit formulae for the derivatives of the two-loop contribution 
to the effective potential involving top, stop, gluon and gluino. In units of a s CfN c /(47r) 3 , where 
Cf = 4/3 and N c = 3 are colour factors, and in terms of the field-dependent quantities defined in 
eqs. (f20|) ~ (j22|) . the 0(a s ) contribution to Vff reads 



AV 



= 2J(m\,m\)-±m\l{m\,m\,G) + 

+ |2 m| I(m| , m| , 0) + 2 L(m| , mf , m?) - 4m t m 5 s 2 g I(m| , mj,m 2 ) 
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, (ci) 



where the functions / , J and L are defined in appendix D. The derivatives of AV that involve only 
the field-dependent stop mixing angle and phase difference (p — fi can be straightforwardly computed 
from eq. (|C1|) . In units of a s Cj?iV c /(47r) 3 , they read 



<9AV 



d 2 AV 



J(m| , m| ) + J(m| , m| ) - J (m| , m| ) 



+2 



m- g m t 
zt dAV 



Tl 2 2 2\ t ( 2 2 2\ 

I(m h ,m- g ,m t ) - I{m l2 ,m~ g ,m t ) 
m g m t 



'26»i 



2 2 2\ tI 2 2 2\ 

m fi ,m g ,m t )- I{m i2 ,m- g ,m t ) 



(C2) 
(C3) 



The explicit expressions for derivatives of AV that involve the quark or squark masses are somewhat 
lengthier. In units of a s CfN c /(4-7r) 3 , they read 



dAV 

dm 2 
ti 



1 if m t \ 2 ( m \ N 

-6 +2mgm t s 2 e t + 4 m t [ 1 - log — ^- ] + 4 m 3 [ 1 - log — f- 



m~ 



m 



ti 



+ 



(5 - c 2 ,^) m\ - s\ Qt m| -Am- g m t s 2dt log 



Q 2 
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/ 9 \ 9 ~ m 

+ (-3 + c 2 " 



4) -|log 2 7# + s2 



+ 



mf 



V log 0^ log ^ + iog g^ iog Qt 

(A + 2 m| mf ) - 2m ^ 20t (mf + mf - m| ) $(m? , mf, mf ) 



2 (mf + to 2 - m? ) -2m- g m t s 2 e, 



m? m? rn 2 " 

T + lo S "TvT 1o § Tvf 

9 



(C4) 



f- -J- (mf + mf - m- g m t s 2 e t ) - s| 0( 



4 

2 

mf 



m 



3 + ^ + — ^ £ 



m 



' mf 



+ 

4irij 

Am? 

ti 

4m| 

Am? 

ti 



( 



2,2 2 
m g + m t — 77l£ 



) 



A -mi 



2(2 2 2 \ 1 
mf + m?) 



A + mf (mf 



m? 



n 

+ (-3 + c 2 et ) log 2 9| + 2(lo g |; 
2 r/ 9 . 9 9\/a « 

A 5 



2 mf 
m~ 



2 

m| 



9 

m~ 










mf + mfj 


log 


mf 

Q 2 






m| 

<5 2 


2 2 Y 
m t -m f J 


log 


m?\ 

















Am? 



mf + mf — m^ ) ( A — 2 mf mf j 



(C5) 



„ / m? \ „ / m 2 

m? (l-log-^f 



'*2 

_mgmt_ 

S20 t 



— m? 
ti 



2 



mr 

1-2 log 



m? m? m? m? 

~ ^2 ("^ t "'/ " "'/ / '■'■■y'ir ' !, :r "'1 



m-gs 2 e t r 2 22 
1 r- 2 - m,- — m 5 — m t 



+ 
+ 



4 

A L 

2 
A 



2 1 2 2 \ 2 2 2\ 

i~ g + m t - m h j $(m fi , m § , m t , 
■,| — mf — mf + 2 mj 



1 1 m § 

L""*i ""* ' i "'" t m 9 S 26»tJ & Q2~ 

2 

22 ( 2 2 2 \ 1 i 

2m g m t - m t m g s 2 e t [m g + m t - m f Jj log ^ 

02/2 2 2 \ m S s 26» ; 
2 m^ lidj — m t — 77i£ J 2 



m t 



(A - 2 mf 



2/2 2 2 

m t (m t — m^ — m^ 



+ f- 2 + m g fl2 * 

W-g S29 t 



m 



1 ti 
log — ± 



mf 



"*t m i 
mf l0g ^ + 1 ° g Q^ l0g ^ 



2 

m t z 



(C6) 



mf 



+ 



1 



Am? 
+2 (: 



L A ( m | - "»1 - 3 mf ) + 2 mf ((mf - m\ ) 2 - mfj 
m? - m|) 3 + 2m? [a + (2ml - mf) (mf + m|)] J $(m?,m?,mf) , (C7) 
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d 2 AV 
dm 2 dm 2 

d 2 AV 



(dm 2 ) 



2\2 



+ 



2 2 

4e t log log 



(C8) 



bm g m\s 20 m 2 
+ 6 log — 2 



2 m\ 



4m| 
4 m? 



A 



"?7 



777* 



^ . rn t g 2 e t / 2 
m ? H m t 



> f 2 

tl 



777, + 777, r + 



ii 



777i 



(m 2 



2 2,2 
m~ g + 777^ 



l?77? 777 2 0™~™2 



A 



2 m § m ( ; s 2 g t 



777 V 



2mf 



(m £ + rrig — 




m ° m i S2 ^ \ log^iiog^ 

2 777^ J m- Q 2 



2 + 



m g 777^ S 2 6»t 



777? 



2 mf 



log — - log — 

g Q 2 g Q 2 



777? 



A 777 2 



"7-5 S2e t 



A 2 

T 



+ ?77 2 ^777- — 2 777 2 + 777^ ^) A + 777^ ^777? — 777 2 + 777^ ^ 

m| (A + (m? + 777 2 ) (2 m 2 - m| )) - (m| - m 2 ) | $(m| , m|, m 2 ) 



777r 
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^2e t 



-S2e t 



(C9) 



a 2 Ay 
9c 2 - am 2 



+ 



2 ?77 4 S 2 8 t 



2 2 / 2 2\ m l 771 

5m h - Am h log — + ( m§ -3m t ) log -| log — 



Q 2 



771? 



+ 777,- 



7777 



777? 



777? 



' os 4 log # +log # log # 



+ 
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t-2 ■ 



S2e t 



— o 1 2 2 1 2 
o — 2 ( 777 = — 777 t + 777 



$(777^, m~,m 2 ) 



S26t 



(CIO) 



where Q is the renormalization scale at which the DR parameters entering the one-loop part of the 
corrections are expressed, the function &(x,y,z) is defined in appendix D, and we have used the 
shortcut A = A(tt7 2 i , m?, ?77 2 ), where the function A(x,y,z) is also defined in appendix D. We recall 
that the derivatives of A^ are computed at the minimum of the effective potential, therefore the 
r.h.s. of eqs. (jC2p (|C10j) is expressed in terms of field-independent parameters (including the mixing 
angle Ot, with —tt/2 < 9t < vr/2). Finally, the derivatives of AV that involve m? can be obtained 



from eqs. (|C4p - (|C7j) by means of the replacements my <-> m^ and s 2 g t — ► —s 2 e t - 
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Appendix D: two-loop functions 



We provide in this appendix the explicit formulae for the two-loop functions appearing in the 0(ata s ) 
corrections to the Higgs mass matrices. Differently from the approach of ref. [18], we choose to 
renormalize the two-loop effective potential before taking its derivatives. As first shown in ref. |36j, 
this is equivalent to using the "minimally subtracted" two-loop functions: 



J(x,y) 



x y 1 - log x 1(1 - log y) , 



(Dl) 



I(x,y,z) 



(x — y — z) log y log z + (y — x — z) log x log z + (z — x — y) log x log y 

A(x,y,z) 



(x + y + z) + 2 [x log x + y log y + z log z 



2z 



(D2) 



L(x,y,z) = J(y,z) - J(x,y) - J(x,z) - (x - y - z)I(x,y,z) . 



(D3) 



In the above formulae, logx stands for log(x/Q 2 ), where Q is the renormalization scale. The functions 
A and $ read, respectively, 



A(x,y,z) = x 2 + y 2 + z 2 - 2 (xy + xz + yz) , 



(D4) 



$ (x,y,z) 



2 logx + logx_ — log u logv — 2 ( Li2(x + ) + Li2(x_) ) + — 



7T" 



(D5) 



where Li2(z) = — Jq dt [log(l — t)/t] is the dilogarithm function and the auxiliary (complex) variables 
are: 

(D6) 



x y 
u = — , v = — 

z z 



A = y (1 - u- u) 2 -4tit), x± = - [1 ± (it - v) - A] 



The definition (|D5P is valid for the case x/z < 1 and y/z < 1. The other branches of $ can be obtained 
using the symmetry properties: 



$ (z, y, z) = * (y, x,z), x $ (x, y, z) = z $ (z, y, x) 



(D7) 



Finally, the following recursive relation for the derivatives of proves very useful for obtaining compact 
analytical results: 



. , . d&(x,y,z) . . . z 

A{x, y, z) = (y + z-x) ®{x, y,z) + - 

ox x 



(y — z) In — V x ( In — + In — 

y V y z 



(D8) 



The derivatives of $ with respect to y and z can be obtained from the above equation with the help 
of the symmetry properties of Eq. (|D7p , 
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